Optical properties of photonic crystal slabs with asymmetrical unit cell 
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Using the unitarity and reciprocity properties of the scattering matrix, we analyse the symmetry 
and resonant optical properties of the photonic crystal slabs (PCS) with complicated unit cell. We 
show that the reflectivity is not changed upon the 180°-rotation of the sample around the normal 
axis, even in PCS with asymmetrical unit cell. Whereas the transmissivity becomes asymmetrical 
if the diffraction or absorption are present. The PCS reflectivity peaks to unity near the quasigu- 
ided mode resonance for normal light incidence in the absence of diffraction, depolarisation, and 
absorptive losses. For the oblique incidence the full reflectivity is reached only in symmetrical PCS. 

PACS numbers: 42.70.Qs, 42.25.Bs 



The physics of photonic crystal slabs (PCS) Q, |1| 
receives mtich attention in recent years because of many 
interesting possibilities they open to control the light in- 
teraction with matter. The tendency is toward the PCS 
sophistication: nanostructured metals or semiconductors 
are included, and the unit cell geometry, using modern 
technology, becomes more complicated. As a result, new 
tools for photonic engineering become available. A well 
known example is the extraordinary optical transmission 
through sub- wavelength hole arrays in metal films 0, 0| . 
Another promising example are polaritonic crystal slabs 
with nanostructured semiconductors 0, im The 
physics behind is the coupling between photonic and dif- 
ferent electronic resonances in such structures. They 
manifest themselves via pronounced resonant Wood's 
anomalies 0, 0, 0| in the optical spectra, due to the 
excitation of qu asig uid ed [3, y, 0, or sur- 

face plasmon P. Il6l ll^ Il8l| or both ^3 modes. On the 
other hand, making complicated unit cells, e.g., charac- 
terized by a lack of 180°-rotational symmetry in the PCS 
plane adds new ways to control the interaction 

with light. 

Thus, the understanding of the symmetry and reso- 
nance properties of the optical response in PCS with 
a complicated unit cell becomes important, also for 
their optical characterization. Measuring reflection from 
asymmetric structure appears to be not a promising 
method for the PCS optical characterization, and a ques- 
tion arises, which optical properties are more sensitive to 
the PCS structure? Meanwhile, notwithstanding a long 
history of the investigations (see, e.g., in [lli. i2L]). start- 
ing actually from the optical gratings, which can be un- 
derstood as one-dimensional (ID) PCS, their properties 
look sometime amazing and even contradictory. 

One example is the nontrivial symmetry properties 
of reflection and transmission from asymmetric PCS. 
Namely, it was demonstrated |l7| experimentally and 
numerically that the reflectivity from 180°-rotationally 
non-invariant (in the PCS plane) metal gratings on a di- 
electric substrate is always symmetric, whereas the trans- 



mission is not. However, the transmission appears to be 
independent on the side of the illumination, whether it is 
from the air or from the substrate (see also in Refs.B|,[i|). 
The authors 17] find this astonishing^ because the cal- 
culated fields distributions inside the PCS appear to be 
very different for the illumination from different sides. 

Another example is the so called anomalous full re- 
flection in zeroth diffraction order in transparent PCS, 
when reflectivity peaks to unity resonantly 'llj. It is es- 
tablished that excitation of surface or quasiguided modes 
in PCS is responsible for this resonant behavior. Such 
resonances, being Fano-type discrete states in the con- 
tinuum background, are characterized by a finite fre- 
quency linewidth due to radiative losses even in trans- 
parent materials. The existence of radiative losses seems 
to exclude a possibility of full reflectivity, however, all 
the models 111 [11 Hi HI H 13 show this effect in 
transparent PCS. Moreover, it appears that very differ- 
ent physical models, with only one common thing, the 
existence of any resonance, predict a qualitatively simi- 
lar behavior of the PCS optical response, see, e.g., dis- 
cussions in [llllll2i.l27t. 

The most general reasons of such behavior are hidden 
in the most general properties of the scattering matrix 
operator of the PCS. These general properties are actu- 
ally well known for many years, but, to the best of our 
knowledge, their consequences for the photonic crystal 
slabs with complicated unit cell have not been yet in- 
vestigated in detail. In this Letter, we recall the most 
general properties of the S-matrix of the transparent ar- 
bitrary PCS and use them to analyse the symmetry and 
resonance properties of the optical response, including 
the conditions of the anomalous full reflectivity. 

Starting from a general scattering matrix operator of 
a planar periodic system made of transparent materials, 
we note that many scattering channels are closed, e.g., 
the backscattering for inclined light incidence, as well 
as all non-Bragg scattering channels. In order to keep all 
nonzero general S-matrix elements, it is convenient to de- 
fine the scattering matrix operator as a unitary 2N x 2N- 
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FIG. 1: Geometry of light reflection/transmission from PCS. 
Solid and dashed arrows correspond to a direct and time re- 
versed process. The main harmonics are shown only. 

dimensional matrix Su, S^^Su = SuSl^ = I, coupling the 
amplitudes of incoming and outgoing harmonics (main 
and diffracted). Here N — 2 + 7Vopon,air + -^open.sub, the 
augend 2 stands for the main harmonic (in air and sub- 
strate), A'opcn.air and A'opcn.sub are the numbers of open 
diffraction orders (into the air and substrate claddings). 
Due to two polarisation states per each harmonic (for ex- 
ample, s andp, or and cr~), there are 2A'' incoming as 
well as outgoing amplitudes (we assume isotropic air and 
substrate claddings). Su = S'u(a;,k||) is a function of the 
real incoming photon frequency lo and in-plane wavevec- 
tor k|| — {kx,ky) — ^ sin??(cosiy9,sin((5), where djCp are 
the azimuthal and polar angles of light incidence. The 
wavevectors of the diffracted photons are 

k||G = (fca:G, fcyc), ^2:(a)G = ^^(a) + (1) 
fc.G - ±^C.%/c2-fc2^-fc2^, (2) 

where G — {Gx,Gy) is the reciprocal 2d PCS lattice, 
and e = 1 or £sub, depending on the PCS side; kzc is 
real for open diffraction channels. Such unitary scat- 
tering matrix Su can be constructed from the infinitely- 
dimensional "large" scattering matrix S accounting for 
the near-field coupling between all the propagating and 
evanescent harmonics and defined, e.g., in Ref. § 
has to be reduced 0| to a "small" 2N x 2A^-dimensional 
matrix S for propagating harmonics only and transferred 
into any energy-flow-orthogonal basis (e.g., of sp or 
polarisations), we do not dwell into details here. 

The symmetry of the system with respect to the time 
reversal in the case of PCS made of non-gyrotropic trans- 
parent materials [compare with Ref. , means that for 
any solution E{u!, k|| ), H{uj, k|| ) of Maxwell equations for 
electric and magnetic fields, i?*(w, — k||), — if*(ci;, — k||) 
remains a valid solution. This property, as well 
known |29j . leads to the reciprocity between the input 
and output channels. In the definition of the unitary 
Su it makes sense to fix the ordering of the channels in 
such a way that the input channels in S'ti(w,k||) trans- 
fer into the output ones in S'„(tj,— k||) and vice versa. 



upon time reversal, see a scheme in Fig^ In the case of 
s,p-polarisations basis, this can be written as 3D| 

|in) = |in(w,k||)) = |out(tj,-k||)*) EE |(?^t). (3) 

Note the important complex conjugates in Eq. Q. Then, 
the reciprocality of the unitary Su matrix means that |3]| 

Suiu;,^)^Sliu;,-^), (4) 

where stands for the matrix transpose. The most gen- 
eral form of Su allowing for the time reversal can be then 
written as 

2N 2N 

Su = ^e*'3^|out)„(in| = ^ e*^^' |out)jj(out|, (5) 

where the input (output) orthogonal bases \va)j (|out)j) 
and the scattering phases Pj, j = 1, . . . 2N are functions 
of (ijJ,k||) which are characteristic for the given PCS. 
The phases are actually even functions of k|| , (3j (w, k|| ) = 
/3j(a;, — k||). Equation ^ means the existence of a spe- 
cial "diagonal" input and output basis sets, which do not 
mix up in the process of scattering. 

It follows from Eq. Q that, with the change of sign 
of k|| in any asymmetric PCS, the reflection with linear 
polarisation conservation is always symmetric, rss(k||) — 
rss(— k||), rpp(k||) = rpp(— k||) (r is the amplitude reflec- 
tion coefficient). Whereas the reflection with depolari- 
sation from p to s (if exists) is asymmetric r<,p(k||) ^ 
rsp(— k||), but equals to that from s to p, rsp(k||) ~ 
rps(— k||). Simultaneously, the transmission (for k||) with 
the conservation of linear polarisation from the air to 
substrate equals always to that in the inverse direction 
(for — k|| and from the substrate to air). 

To the contrary, in the case of circular polarisations, 
because now the time reversal switches from to u^, 
the reflection with change of polarisation is symmetri- 
cal, r^ (k||) — r_| (— k||). Simultaneously, r++(k||) = 

r (— k||), but the reflection with conservation of polari- 
sation is asymmetrical, ^^^(kii) ^ r_|_+(— k||). Addition- 
ally, 7'+_(k||) ^r_+(-k||). 

These symmetry properties of the reflection (and 
transmission in reverse direction) are independent of how 
many diffraction channels are open. They are even more 
general and hold for PCS with absorptive materials, sim- 
ply because any absorption channel can be included into 
a more general unitary S-matrix of the full system as an 
additional scattering channel. The symmetry of trans- 
mission in reverse direction is just what was found in 
Refs. There is no contradiction here with the 

different field distributions, because the reciprocity exists 
between the time-reversal channels only. Whereas the full 
solution is not reciprocal, due to different reflected waves 
in the cases of incidence from the air and the substrate. 

As to the symmetry of transmission for illumination 
from the same side of the PCS, in case of asymmetric 
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FIG. 2: Calculated angular dependences of transmissivity T — \tf (a,d,g), reflectivity R = \rf (c,f,i), and total diffraction 
D = 1 — T — R (b,e,li) spectra of the s-polarised light in the asymmetric (a-f) and symmetric (g-i) ID PCS with cross-sections 
and parameters shown in the inserts. For panels (d-f) and (g-i) the angles of incidence i} are indicated in the legends in panels (f) 
and (i), respectively. Open triangles in panels (g,i) show the resonance model Eas. lll3ll4ll . the resonance parameters indicated 
in panel (h). The plane of s-polarised light incidence is perpendicular to grooves {^p = 0). The dips (peaks) in transmissivity 
(reflectivity) are quasiguided modes, a less pronounced wrinkles (or cusps) are the diffractive Wood's anomalies. Note that the 
reflectivity is symmetrical, the diffraction is asymmetrical, and the transmissivity is symmetrical if only the diffraction is zero. 



PCS it holds at frequencies when aU diffraction channels 
are closed, and only in transparent PCS, because the k|| 
and — k|| transmission processes from the same side of 
the PCS are not reciprocal, see in Fig^ 

We give an example of such asymmetric optical re- 
sponse in Fig. 121 a-f ) of a model transparent asymmetric 
ID PCS, see the schematic crossection in panel (e). The 
incoming light is s-polarised (electric field parallel to the 
PCS grooves). The calculations are done via §-matrix 
method from Ref. using 21 harmonics in Eq. 0J with 
accuracy better than 10~^. Note that so far the diffrac- 
tion is absent (for photon energy below w 1.65 - 1.75 eV, 
depending on the angle of incidence) the transmissivity 
is symmetrical with the sign of k||. Note also that, in 
agreement with the discussion above, the reflectivity is 
always symmetric, see in panels (c,f). 

The transmissivity dips in Fig. El are Wood's anoma- 
lies 0, 0, 0| due to quasiguided modes in PCS. Near 
such resonances the reflectivity can become full in the ze- 
roth diffraction order. Numerical examples of such a full 
reflectivity can be seen in Figs.[2t,f for normal incidence 
in asymmetric ID PCS, and for any z9 in a symmetric 
PCS, Fig. 121. The full reflectivity is achievable unless 
the first diffraction channel opens, compare with panels 
(b,e,h). 

These properties of the resonant reflectivity in trans- 
parent PCS can be deduced from the general properties 
of the unitary scattering matrix Su- Let detS~"'^(ci;,k||) 
has a zero at wo(k||) = rio(k||) — i7o(k||) in the lower half 



of complex energy plane, corresponding to the guasigu- 
ided photonic mode 



-i(c.o,k||)|0(k||)) =0, 



(6) 



where |0) is the resonant output eigenvector in the 
"large" basis. Here we analyse the simplest case of a 
nondegenerate quasiguided mode. 

Su has to be unitary for real lo, and it imposes sig- 
nificant restrictions on its possible form. In fact, near 
the resonance one of the scattering phases (the reso- 
nant phase) is a quickly changing function of energy, and 
Eq.llSJ can be rewritten in a form 

5„(t^,k||) = ^e'ft |o)„-(o| +77(60, k||) |o)ii(o|, (7) 

w-wo(k||) w-wo(k||) 

where |o)i = |o(k||))i, i.e., the resonant "small" outgo- 
ing vector in the orthogonal basis corresponding to the 
"large" resonance vector |0), |6)i = |o(— k||)*)i, and 
\o)j^i are the output basis vectors in the subspace or- 
thogonal to |o)i. To ensure 77(0;, ky) = 77(0;,— k||) (the 
reciprocity of the reflection near the resonance), we have 
to set (k||) = ti;o(^k||). However the resonant vec- 
tors |o(k||))i and [o(— k||))i can be different if there is 
no additional symmetry of the structure transforming k|| 
to — k||. Equation Q is the most general form of the 
Breit-Wigner formula (see, e.g., in textbook [2^]) for the 
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resonant optical response in a transparent PCS. Except 
r], all quantities in Eq. ^ are slow functions of uj; ne- 
glecting this UJ dependence gives a very good resonant 
approximation for in the vicinity of fig |32l |. 

In the zeroth diffraction order Eq. ||7J) can be simplified 
further. In this case Su is a 4 x 4 matrix. However, in the 
case of ID PCS, and for linearly polarised light with plane 
of incidence perpendicular to the grating grooves (ip — 
0), the s and p polarisations are decoupled. Then, Su 
becomes effectively a 2 x 2 matrix. By a proper selection 
of the outgoing harmonics phases, the corresponding 2D 
vector |o)i can be made real, and 

i°H=(::4)'i°>-(-iy- (') 

The most general form of the (2 x 2)-dimensional unitary 
scattering matrix in this polarisation becomes 

5„(c^,k||) = (cos|,-sin|)e^'' 

+ (^'^"^^) (sini, cosDrKc.,k||), (10) 

where /3,^, and ^ = ^(— k||) are slowly changing with 
UJ, k|| parameters of the system near the resonance. For 
the coefficients of reflection r = {Su)ii and transmission 
t = {Su)2i we have, respectively 

r{uj, k||) — e^^ cos^ cosf + r]{uj, k||) sinf sinf . (11) 
t{u), k|| ) ~ — e'*^ sin f cos ^ + r]{uj, k|| ) cos f sin f . (12) 

In the case of normal incidence (k|| — 0) for any PCS, 
we have obviously that ^ = C and, as a result, 

r(w,k||) = e*'^cos^C + ^('^,k||)sin2^, (13) 
t(w,k||) = [r/(tj,k||) -e*'3]cosCsin^. (14) 

Obviously t = and \r\ = 1 at the energy when r]{uj) = 
e'^ . Because /3 is a slow function of energy, and rj makes 
a circumnutation from ■i]{uj ^ f^o — 7o) ~ — 1 through 
77(f2o) = 1 to r]{uj ^ rig —70) ~ ~1: this condition always 
matches at some energy near the resonance, provided e^^ 
is not too close to -1. 

The quantity sin ^ cos ^ = 1/2 sin 2^ reaches its maxi- 
mum 1/2 at ^ = 7r/4. Thus, the only possibility to reach 
a full transmission |i| = 1 is to have ^ = 7r/4 and, as seen 
from Eq. H14|l . ri{u)) = —e^^. For example, the full trans- 
missivity is reached exactly at the resonance, w = fioj in 
the only case if /3 = tt and ^ = 7r/4 simultaneously. This 
is the case of a well known Fabry-Perot resonator with 
symmetrical lossless mirrors. 

In asymmetric PCS, at oblique incidence, as can be 
understood from the inspection of Ea. (|ll|l . the reflec- 
tivity still peaks (or anti-peaks) near the resonance en- 
ergy. However, now the maximum value of |r| is obvi- 
ously less than 1, because £, ^ ^, i.e., the resonant vectors 
o(k||))i 7^ |o(— k||))i for the asymmetric PCS. 



In case if the PCS has an additional symmetry trans- 
forming the input channels with k|| to that with — k|| , we 
have £, = ^. Then Eqs. (|13I14(I hold for any angle As a 
consequence, the reflection can be full near the resonance 
for oblique incidence, not only for k|| = as in the asym- 
metrical gratings. This additional symmetry can be, e.g., 
a vertical mirror plane (symmetrical PCS on a substrate, 
see insert in Fig. The numerical examples shown in 
Fig. EJi-f fully agree with this analysis. We would like to 
add here that the property of full reflectivity is obviously 
destroyed by any losses, including absorption, defect and 
PCS edge scattering for finite sample size. 

The resonant approximation neglecting slow uj depen- 
dences in Eqs. 11HI - (|14|I gives a very good description of 
the optical response near the resonance, see an example 
in Fig. E^-i. The parameters flQ,jo,^,^, and /3 for any 
PCS can be calculated from the full scattering matrix 
§ ,32.] . This resonance approximation gives, e.g., a rigor- 
ous explanation of an intuitive model of the interference 
between direct and indirect pathways [l3 |. with minimal 
number of directly calculated parameters (5 for general 
PCS, and 4 for normal incidence or symmetric PCS), and 
it is applicable for arbitrary PCS. 

To conclude, using the unitarity and reciprocity prop- 
erties of the scattering matrix, we analyze theoretically 
the nontrivial symmetry properties and near resonance 
behavior of the optical response in photonic crystal slabs 
(PCS) with asymmetric unit cell. As a direct conse- 
quence of the reciprocity, the reflection with conserva- 
tion of linear polarisation is always symmetrical, whereas 
that with depolarisation, if exists, is asymmetrical. For 
the circularly polarised light the opposite rule holds. As 
a direct consequence of the unitarity, the PCS reflectiv- 
ity peaks to unity near the quasiguided mode resonance 
for normal light incidence in the absence of diffraction, 
depolarisation, and losses. For the oblique incidence the 
full reflectivity in zero diffraction order is reached only in 
symmetrical PCS; in asymmetrical PCS the reflectivity 
maximum decreases with the angle of incidence increase. 
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